J/iJj elliptic flow in relativistic heavy ion collisions 

Yunpeng Liu a , Nu Xu b and Pengfei Zhuang a 

a Physics Department, Tsinghua University, Beijing 100084, China 

b Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 
94720, USA 

The J/ij} elliptic flow in high energy nuclear collisions is calculated in a transport model. 
While the flow is very small at SPS and RHIC energies, it is strongly enhanced at LHC 
energy due to the dominance of the regeneration mechanism. 

1. introduction 

The J/ip production is commonly accepted as a signature of the quark-gluon plasma 
formed in high energy nuclear collisions. The J/ip elliptic flow which is sensitive to the 
early space-time evolution of the colliding system carries dynamic information on the new 
state of matter, if it is formed in the early stage. Since the recently observed J/ip elliptic 
flow at RHIC energy [[TJ is with large uncertainty and can not differentiate from various 
J I ip production mechanisms, it needs more precise experimental and theoretical studies. 

The charm quark number in the partonic phase increases with the colliding nuclear 
energy. When the energy is high enough, the coalescence of those uncorrelated charm 
quarks becomes a significant source for J/ip production. The J/ip regeneration and its 
competition with the primordial J / ip production have been discussed in different models [ 
[21 El 0]. To calculate the fractions of the initially produced and regenerated J ftps for 
the nuclear modification factor and transverse momentum distribution, we set up the 
charmonium transport equations [[5] with both gain (regeneration) and lose (suppression) 
terms in the partonic phase. In this paper, we calculate the J/ip elliptic flow in nuclear 
collisions at SPS, RHIC and LHC energies in the frame of the transport model. 

2. transport model 

Taking the fact that charmonia are heavy particles, the distribution function fe(p t , x t , r|b) 
where \l/ stands for J/tp,x c and ip' in transverse phase space at mid-rapidity can be de- 
scribed by a Boltzmann equation [[5] 

dfo/dr + • V/* = -amfm + fa, (1) 

where p f ,x t ,r and b are respectively the charmonium transverse momentum, transverse 
coordinate, invariant time and impact parameter, v# = Pt/Ey is the charmonium trans- 
verse velocity, a and (3 are dissociation and regeneration rate calculated from the gluon 
dissociation process \l/ + g — > c + c and the reverse process. From the large single electron 
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elliptic flow measured at RHIC [[6], charm quarks are supposed to be kinetically ther- 
malized in the hot medium. The temperature, baryon chemical potential and medium 
velocity appeared in the gluon and charm quark thermal distributions are determined 
from the 2+1 dimensional ideal hydrodynamics. Solving the coupled transport equation 
for J/ip and the hydrodynamic equations for the space-time evolution of the partonic 
phase, one can obtain the J/ip distribution function fj/^(pt , x 4 , r|b) and in turn the 
nuclear modification factor Raa and the averaged transverse momentum square [|3|. 

3. J/ip elliptic flow 

The initially produced charmonia from nucleon-nucleon hard processes carry high mo- 
mentum, and the gluon multi-scattering effect before the two gluons fuse into a charmo- 
nium leads to a further transverse momentum broadening [[?]. Therefore, the high p t 
J/ips in the final state are dominated by the initial production. On the other hand, the 
thermalized charm quarks in the hot medium satisfy the thermal distribution, and the 
regenerated charmonia are with low momentum and thus dominate the J/ip production 
at low p t region. The left panel of Fig{T] shows the fractions of initially produced and 
regenerated J/ips as functions of pt for central Pb+Pb collisions at LHC energy. It is 
easy to see that the low pt region is governed by regeneration and the high p t region is 
characterized by initial production. 

The elliptic flow v 2 describes the asymmetric degree of the particle momentum in trans- 
verse space, v 2 (pt) = (pl — V^l/Vt- At SPS energy, there are only few charm quarks 
produced at the initial stage of the collisions. Therefore, the regeneration can be safely 
neglected. In this case, the initial production controls the system, and the non-thermalized 
J/ips can not feel the collective flow of the thermalized medium. The only source for the 
J/ip elliptic flow is the leakage effect described by the free streaming term on the left 
hand side of the transport equation (jTJ): J/ips with high p t are easier to escape the hot 
medium in the direction where the fireball is thinner. Therefore, the final state J/ips are 
no longer isotropic and leads to finite value of v<i. However, due to the expansion of the 
fireball, the asymmetry in geometry decreases with time, and the elliptic flow due to such 
geometry configuration is very small in comparison with the collective effect. We show 
on the right panel of Fig{T] the elliptic flow in heavy ion collisions at impact parameter 
b = 7.8 fm. At SPS, the elliptic flow is extremely small due to the lack of regeneration. 

The fraction of regeneration for J/ip increases with centrality monotonously. At RHIC, 
the regeneration becomes important and its contribution to the total J/ips is even around 
50% in central Au+Au collisions [@J IS]- However, for semi-central collisions where the 
elliptic flow reaches the maximum, the system is still controlled by the initial production 
and the J/ip elliptic flow is still very weak, as showed in the right panel of Fig{TJ 

The situation at LHC is very different. The formed fireball is much larger, hotter and 
longer lived, almost all the initially produced J/ips at low and intermediate momentum 
are eaten up by the medium, and the J/ip production with p t < 4 GeV is dominated by 
the regeneration, see the left panel of Fig{TJ In this case the J/ips from the recombination 
of thermalized charm quarks carry large elliptic flow, and the maximum value at about 
Pt ~ 4 GeV reaches 10%. We have assumed charm quark thermalization at any p t , the 
elliptic flow from the regeneration thus increases monotonously with p t . While this is 
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not true at high p t , it does not change the full v% remarkably, since the high p t region 
is dominated by the initial production. Note that the shape of the J /if) v 2 (pt) is quite 
different from what observed for light quark hadrons [[5] which is saturated at high p t . 

In summary, the J/if> elliptic flow at LHC is expected to be much larger than that 
at SPS and RHIC, due to the dominance of the regeneration mechanism and the full 
thermalization of charm quarks. 




p t (GeV/c) P t ( QeV/c ) 

Figure 1. The J/if> production fractions (left panel) and elliptic flow (right panel) as 
functions of transverse momentum in Pb+Pb collisions at impact parameter b = 7.8 fm at 
LHC energy. The calculations with only initial production, only regeneration (multiplied 
by a factor of 0.2) and both are indicated respectively by dashed, dotted and solid lines. 
The full elliptic flow at SPS and RHIC are also showed as a comparison. 
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